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ABSTRACT 

Aims. AGN-driven outflows are invoked by galaxy evolutionary models to quench star formation and to explain the origin of the 
relations observed locally between super massive black holes and their host galaxies. This work aims to detect the presence of 
extended ionised outflows in luminous quasars where we expect the maximum activity both in star formation and in black hole 
accretion. Currently, there are only a few studies based on spatially resolved observations of outflows at high redshift, z > 2. 

Methods. We analyse a sample of six luminous (L > 10 47 erg/s) quasars at z ~ 2.4, observed in H-band using the near-IR integral 
field spectrometer SINFONI at VLT. We perform a kinematic analysis of the [Om] emission line at A = 5007A. 

Results. We detect fast, spatially extended outflows in five out of six targets. [Om]/I5007 has a complex gas kinematic, with blue- 
shifted velocities of a few hundreds of km/s and line widths up to 1500 km/s. Using the spectroastrometric method we infer size of 
the ionised outflows of up to ~2 kpc. The properties of the ionised outflows, mass outflow rate, momentum rate and kinetic power, 
are correlated with the AGN luminosity. The increase in outflow rate with increasing AGN luminosity is consistent with the idea that 
a luminous AGN pushes away the surrounding gas through fast outflows driven by radiation pressure, which depends on the emitted 
luminosity. 

Conclusions. We derive mass outflow rates of about 6-700 M© yr _1 for our sample, which are lower than those observed in molecular 
outflows. Indeed physical properties of ionised outflows show dependences on AGN luminosity which are similar to those of molecular 
outflows but indicating that the mass of ionised gas is smaller than that of the molecular one. Alternatively, this discrepancy between 
ionised and molecular outflows could be explained with different acceleration mechanisms. 

Key words, galaxies: active - galaxies: evolution - quasars: emission lines - techniques: imaging spectroscopy 


1. Introduction 


Feedback mechanisms from quasars (QSQ s) are considered to 
be crucial for galaxy evolution (see King 2010| |Fabian|[2012| 


and |King & Pounds||2015| ). During the bright active phase, Ac 


tive Galactic Nuclei (AGN) are believed to drive energetic out¬ 
flows that expel gas at large scales from their host galaxies 


(e.g Granato et al. 2004, [Di Matteo et al. 

2005 Menci et al. 

2008 King 2010[]Zubovas & King 2012 Fal 

bian 2012, Faucher- 


2014||Costa et al.|20i4||Costa et al.|2015| ), hence removing the 


* based on Observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, RID: 
086.B-0579(A) 


supply of cold gas required for star formation (SF) activity. Ac¬ 
cording to some of these models, the black hole (BH) achieves 
a critical mass during the later stages of a merging event and the 
energy output of the associated AGN is so large that radiation 
pressure drives a fast outflow in the nuclear region which sweeps 
away the gas in the host galaxy. When the feedback phase is 
over, the stellar population in the host galaxy continues to grow 
mainly through minor or major mergers with other galaxies, (e.g. 
|Baldry et al.l2004[|Perez-Gonzalez et al.|2008[ ). At the same time 
the lack of gas around the nucleus slows down the growth of the 
BH. Models invoke feedback mechanisms to explain the origin 
of the correlation between the mass of supermassive BHs and the 
mass and velocity dispersion of host galaxy bulges observed in 
the local Universe (e.g. [Magorrian et al.|19~98j |Marconi & Hunt 
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|2003| |Ferrarese & Ford|[2005[ |Kormendy & Ho||2013| ). Essen¬ 
tially, the observed correlations can be explained as the result 
of the balance between the outward radiation force generated by 
the AGN and the inward gravitional force of the host galaxy (e.g. 
|King|2010[|Labian|2012[ ). 

Quasar-driven outflo ws extending to kpc-scales have been 
resolved both locally (e.g Feruglio et al. 2010} Rupke & Veilleux 


2011; Cicone et al.j|20 12 "[Rupke & Veilleux 2013; Rodriguez 

Zaurin et al.||2013[ |Leruglio et al.||2013b|a| |Cicone et al.||2014| 


Rodriguez Zaurm et al. 2014; Aalto et al. 2015; Leruglio et al. 

2015 

and at high redshift (e.g. Alexander et al. 2010; Nesvadba 

et al. 

2011; Harrison et al. 2012!; 

Maiolino et al. 2012 

Cano- 

Diaz et al. 2012; Harrison et al. 20; 

L4; Cresci et al. 2015; 

Cicone 


et al.|2015p in ionised, atomic and molecular gas. Despite big ad 
vances in data quality and analysis in the past decade, the main 
properties of these energetic outflows remain largely unknown. 
The exact outflow morphologies and driving mechanisms are 
still poorly known: it is debated whether their morphology is 
conical or shell-like, and it is also unclear what is the physical 
process responsible for the coupling of the energy/momentum 
released by the central AGN with the galaxy interstellar medium 
(e.g. inner winds and shocks, radiation pressure on dust). Ac¬ 
cording to one of such scenarios, the fast wind, accelerated close 
to the BH by radiation pressure, shocks the ISM of the host 
galaxy creating a bubble of hot gas which expands at large ve¬ 
locities (~ 1000 km/s). If the post shock material does not cool 
efficiently, energy is conserved, the bubble expands adiabatically 
and the outflow is energy-driven. On the other hand, if the post 
shock material cools efficiently, due for example to Compton 
cooling by AGN photons, only momentum is conserved and the 
outflow is momentum-driven. Since the efficiency of the cooling 
process, dominated by inverse Compton scattering, drops with 
increasing shock radius, there is a critical distance (~ 100 pc - 
1 kpc) beyond which the extended outflows can only be energy- 
driven ( [Zubovas & King||2012| ). Indeed, recent observati ons of 


kpc-scale outflows in local AGNs support this scenario (Cicone 
|et al.|2014||Feruglio et al.|20T5)|Tombesi et al.|2015| ). 


Moreover it is not clear whether molecular and ionised out¬ 
flows are accelerated by the same mechanism, whether they have 
the same spatial distribution or whether they occur on similar 
timescales. So far, there are only a few observations of both 
molecular and ionised AGN-driven outflows in the same galaxy: 
SDSS J1356+1026 is an example of an obscured QSO where 
molecular and ionised outflows have different properties, i.e. 
outflow rates, velocities, radii, morphologies and time scales 
( [Greene et al.|2012[|Sun et al.|2014| ). 

The forbidden emission line doublet [Oim] at A = 
5007,4959A is a good tracer of ionised outflows on large scales 
since it can not be produced at high densities and so it cannot 
trace the sub-parsec scales of the Broad Line Region (BLR). In 
the presence of outflows, the spectral profile of the [Oiii]/15007 
emission line can be highly asymmetric, with a broad, blue- 
shifted wing that is rarely observed in star-forming regions. Lor 
this reason, the [Om]/I5007 emission line has been used to iden¬ 
tify outflowing ionised gas in low-redshift and high-z AGNs. Re¬ 
cent integral field observations of the [Om]/I5007 emission have 
provided quantitative measurements of the outflowing gasprop- 
erties in AGN at low redshifts z < 0.5 (e.g. Greene et al. |2011| 
Harr ison et al.|2014| ). At higher redshift, z ~ 2, [Alexander et ah] 


( 2010| ) found evidence for galactic-scale ionised outflows by 
mapping the [Oiii]/15007 emission line in ULIRGs (Ultra Lumi¬ 
nous IR Galaxies) hosting type 1 AGNs. [Cano-Diaz et al.| ( |2012] ) 
observed a ionised outflow, extended up to ~8 kpc, in a QSO 
at z ~ 2.4. Linally, |Brusa et al.] ( |2015| ) and |Perna et al.| ( [2015] ) 
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detected broad [Oiii]/15007 blue wings in a sample of obscured 
AGN revealing the presence of outflows extended over several 
kiloparsecs. In one case of the latter sample, the ionised outflow, 
extended up to ~ 12 kpc, is anti-correlated with the presence of 
star formation in the host galaxy (Cresci et al. 2015) ), similarly 
to what found by |Cano-Diaz et al.| ( |2012| ). 

In this paper, we present a kinematical analysis of the 
[Oiii]/15007 emission line observed in six high luminous (Lb 0 i > 
10 47 erg/s) quasars at z ~ 2.4. The line profiles and the velocity 
maps obtained by the kinematical analysis show the presence of 
ionised outflows extended on scales larger than 2-3 kpc from the 
nucleus. This is one of the first observations of extended ionised 
outflows in QSO or type 1 AGN; indeed most AGN-driven out¬ 
flows have been observed so far in type 2 AGN. We present a new 
method to measure outflow properties from ILU (Integral-Lield 
Unit) data when the source is only marginally spatially resolved. 
We show that the typical signature of outflows, i.e. asymmetric 
line profiles and blue-shifted components in the velocity maps, 
can be described by two point sources (the central AGN source 
and the outflowing material) separated by a distance of a few 
kpc. We find that the ionised gas likely traces only a fraction of 
the total outflowing gas, unless ionised outflows have a different 
origin than the molecular ones. 

The paper is organised as follows: in Section 2 we present the 
sample selection and properties, in Section 3 we show the data 
analysis and spectral fitting. The results of our data analysis are 
presented in Section 4. Linally, in Section 5 we discuss the main 
results, i.e. the nature of ionised outflows and the comparison 
with other observations. A Ho = 67.3 km s -1 Mpc -1 , Qm = 
0.315, = 0.685 cosmology is adopted throughout this work 

( [Planck Collaboration et al.|2014| ). 


2. Sample selection 

We selected six QSOs at z ~2.4 with large [Oiii]/15007 equivalent 
widths (> 10A in the rest frame) and bright H-band magnitude 
(< 16.5 mag) from the sample of Shemmer et al. ( |2004| ), Netzer 
|et al.| ( |2004] ) and |Marziani et al.| ( 2009| ). The luminosities of the 
objects in our sample are in the range L = 10 47 -10 48 erg/s, mak¬ 
ing them the highest luminosity sources where outflows have 
been spatially mapped. These characteristics are chosen to max¬ 
imise our chances to detect signatures of feedback on the host 
galaxy. As explained in Section[l] the cosmic epoch correspond¬ 
ing to z - 2 is crucial for the growth of the most massive galax¬ 
ies and black holes, and shows the maximum activity both in star 
formation and in BH accretion (e.g., |Madau & Dickinson|2014| . 
The selected objects, at the high end of the QSO luminosity func¬ 
tion, are those where we expect to detect feedback in action, at 
least according to current galaxy evolutionary models (e.g. |Hop-| 
[kins et al.|2006| . Moreover, the targets were selected depending 
on their particular [Om] profile. Indeed, the large [Om]/15007 
equivalent width should allow us to easily recover the kinemati¬ 
cal maps of the ionised gas, a challenging task in luminous QSOs 
where narrow emission lines are usually weak (e.g Netzer et al. 
|2004| ). Linally, the bright ionised line profile can be easily de- 
blended from the broad H ft and Lell emission, which are associ¬ 
ated to the BLR. Indeed, based on near-IR spectra available from 
the literature ( Shemmer et al.[2004[ [Netzer et a l.|2004 ; Marziani 
|et al.|2009) , we also chose our objects such that the contamina¬ 
tion from broad Lell emission is expected to be small. 

Basic properties of the observed objects are given in Table [T] 
One of our sources, 2QZJ0028, is already studied by Cano-Diaz 
et al.| ( |2012| ), and was reobserved by using much longer integra¬ 


tion time, confirming their results. 
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3. Observations and data reduction 


The targets, whose properties are described in Table [T] were 
observed in February and September 2012 using the near-IR 
integral field spectrometer SINFONI (Spectrograph for INte- 
gral Field Observations in the Near Infrared) at the VLT (Very 
Large Telescope). The observations were obtained in H-band 
(A ~ 1.45 - 1.85 p m, where [Oiii]/15007, [Oiii]4959 and H J3 are 
detected), in seeing limited mode with 0.250" spatial scale and 
medium spectral resolution of R = 3000. In all observations, the 
object was moved across the 8 "x 8 " field of view in order to per¬ 
form the sky subtraction. A standard star for telluric correction 
and flux calibration was observed shortly after or before the on- 
source exposures. The total on-source integration time is over 3h 
per target. The airmass are different for each target and spanning 
a range between ~ 1.0 and ~ 1.4. 

After removing cosmic rays from the raw data using the 
L.A. Cosmic procedure (van Dokku m|2001| ), we used the ESO- 


SINFONI pipeline to reduce the data. The final data cubes pro¬ 
duced by the ESO-SINFONI pipeline have a spatial scale of 
0.125"x0.125" and a field of view of about 8 "x 8 ". The esti¬ 
mated angular resolution is ~ 0.4-0.6", based on a 2D-Gauusian 
fitting of the flux map of the spatially unresolved broad H {3 line 
(see Section [471]). 


4. Data analysis 

Figure [T] shows the H-band spectra extracted from the nuclear 
region of each QSO from an aperture of 0.25"x0.25". The spec¬ 
tra clearly show the broad Hj3 and the strong emission line dou¬ 
blet [Oiii]/L15007,4959. In addition to these components, weak 
Fell emission lines are also visible in two out of the six QSOs 
(LBQS0109 and HE0251). The asymmetric [Oiii]/15007 profiles 
suggest the presence of ionised outflows in most of the targets 
(LBQS0109, 2QZJ0028, HB8905, HB8903 and HE0251). In¬ 
deed the presence of a prominent blue wing is the typical sig¬ 
nature of high velocity gas moving toward the line of sight. The 
H (3 emission line shows asymmetric velocity profiles as well as, 
in some cases, two distinct emission peaks (i.e HB8903). Indeed, 
in all spectra the H ft line is a sum of two different components: 
a very broad (FWHM > 2000 km/s) one which is associated to 
the BLR emission and a fainter, narrower ( FWHM~ 500-1200 
km/s) one likely associated to NLR emission. 

4.1. Fitting procedure 

In order to understand the dynamic and the main properties of 
the ionised outflows, we performed a kinematical analysis on 
the forbidden lines. 

At first, we extracted the spectrum of each QSO from a nu¬ 
clear region of 0.25"x0.25" where the signal-to-noise (S/N) is 
highest. We fitted simultaneously the continuum, H ft and [Om] 
emission lines by minimizing x 2 using the IDL routine MPFIT 
( |Markwardt|2009[ ). 

The [Om]T5007 line profile is very complex in these objects 
and a single Gaussian function is not able to reproduce the asym¬ 
metric velocity profiles, so we fitted the emission-line doublet 
using either multiple Gaussian components or a broken power- 
law convolved with a Gaussian distribution. Initially, we fitted 
the emission line using two Gaussian components, one broad 
(FWHM > 1000 km/s) and one narrow (FWHM < 1000 km/s), 
and, when necessary, added a third component to minimise the 
X 2 value. Since in a few cases the [Oiii]/15007 emission line ex¬ 
hibits a very asymmetric line profile with blue wings extended 


to |v| >1000 km/s, we replaced the multi-Gaussian components 
with a broken power-law profile convolved with a Gaussian dis¬ 
tribution. We did not add additional Gaussian components in or¬ 
der to avoid an unnecessary large number of free fitting parame¬ 
ters. We do not attribute any physical meaning to the individual 
fitting component and we measure gas kinematics by analysing 
the total line profile. Therefore, using either multiple Gaussian 
components or a broken-power law does not change the results 
of this work, provided that the quality of the fit is similar. The 
two emission lines of the [Om] doublet, originating from the 
same upper level, were fitted imposing the same central velocity 
and velocity dispersion, with the intensity ratio 7(5007)//(4959) 
fixed at ~3. 

H/? is well described by two components, a broad (BLR) and 
a narrow one (NLR). We used a broken power-law profile for the 
very broad component (FWHM > 2000 km/s) since it usually 
provides a good fit to QSO broad emission lines (Nagao et al. 


12006| . The model used for the H p profile of 2QZJ0028 shows a 
broad red wing that is likely due to continuum or Fell emission, 
so we do not attribute any physical meaning to this red-shifted 
emission. A similar profile is visible in Figure 1 of |Cano- Dfaz 


|et al.| ( [20T2] ). Since the narrower Up is weak and sometimes only 
marginally detected, it is not possible to reliably constrain its 
profile and kinematics and therefore we assume that it has the 
same average velocity and velocity profile as [Oiii]/15007. 

Finally we used a power-law for the continuum and, for those 
spectra showing Fell emission, we used the Fell template from 
Tsuzuki et al.| ( |2006] ). These best-fitting solutions are used as first 
guesses in the pixel-to-pixel fitting of the kinematical analysis 
described in the following. 

Before performing a spatially resolved kinematical analy¬ 
sis, we tested whether the kinematics of the ionised gas is spa¬ 


tially resolved. Following the method described in Carniani et al. 


( |2013| ), we analysed the flux residual maps obtained with a pixel- 
by-pixel kinematical fitting with the components just described, 
after assuming that the targets are spatially unresolved. In this 
case, the H-band spectrum is expected to be the same in any 
spaxel apart for a different normalisation. The residual maps 
(Figure [2]) at the wavelength range of [Oiii]/15007 indicates that 
the forbidden line emission is kinematically resolved in all but 
one of the quasars, i.e HE0251. In fact, if the emission was unre¬ 
solved, we would expect to observe both a spectral profile and a 
residual map consistent with the noise, as in the case of HE0251. 
The flat residual at the wavelength range of H P indicates that the 
broad component of H P is spatially unresolved, consistently with 
a H p origin in the BLR. 

After testing that the [Om]/i5007 line emission is kinemat¬ 
ically resolved in five of the six sampled QSOs, we performed 
a pixel-by-pixel fitting on the data cubes allowing the spectral 
components to vary except for the broad Up which is not spa¬ 
tially resolved. Figure [3] shows the kinematical properties of the 
ionised gas as obtained from the fit of the [Oiii]T5007 line: 


o zeroth moment map (i.e., line flux map; first panel); 
o first moment map (i.e., median velocity map; second panel); 
o blueshifted velocity map, vio, the velocity at the 10 th per¬ 
centile of the overall emission-line profile (third panel); fitted 
in each spatial pixel; 

o line width map , W$ o, the velocity width of the line that con¬ 
tains 80% of the emission-line flux (fourth panel); this is de¬ 
fined as Ws o = V 90 - vio, where vio and V 90 are the veloci¬ 
ties at 10th and 90th percentiles, respectively. For a Gaussian 
profile, Wg 0 is approximately the full width at half maximum 
(FWHM). 
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Table 1 . Properties of our quasar sample 


ID 

Target Name 

[pm] 

Arrest [A] 

Redshift 

AL 5100 [10 46 erg/s] 

LBQS0109 

LBQS 0109+0213 

1.68 + 0.04 

28.5 + 0.2 

2.35 + 0.08 

4.9 ± 0.8 

2QZJ0028 

2QZ J0028-2830 (a) 

1.70 + 0.04 

27.9 + 0.2 

2.40 + 0.09 

3.1 ±0.6 

HB8905 

HB89 0504+030 

1.75 + 0.05 

22.7 + 0.2 

2.48 + 0.09 

2.9 ± 0.6 

HE0109 

HE 0109-3518 

1.706 + 0.003 

17.342 + 0.011 

2.407 + 0.007 

7.4 ± 1.5 

HB8903 

HB89 0329-385 

1.720 + 0.013 

19.31+0.06 

2.44 + 0.03 

5 ± 2 

HE0251 

HE 0251-5550 

1.68 + 0.03 

31.81+0.17 

2.35 + 0.05 

6.8 ± 1.4 


(a): this is the same target of Cano-Diaz et al. ( 2012) which has been re-observed. 
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Fig. 1 . Upper panel. The spectra of the six QSOs if our sample. Each spectrum is extracted from a nuclear region of 0.25" x025" (2x2 pixel). 
The different components in the fit for each line (H/A[Oiii] and Fell) are shown in green and the red line is the total fit. The shaded yellow regions 
indicate the zone affected by strong sky line residuals which are excluded from the fit. Lower panel fit residuals, obtained as a difference between 
observed and model spectra. 


The observed kinematical maps are the results of the convolu¬ 
tion of the intrinsic ones with the PSF of the observations. In 
those cases where the [Om]T5007 emission line has been fitted 
with multiple Gaussians, line profiles are based on the sum of 
all Gaussian components. The maps were obtained by selecting 


only those spatial pixels with a S/N equal to or higher than 2. We 
defined the S/N as the ratio between the peak of the [Oiii]/15007 
line and the rms of the residuals. Zero-velocities correspond to 
the peak position either of the narrow (FWHM < 1000 km/s) and 
strongest Gaussian component or of the broken-power law one 
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Fig. 2. Residuals from the pixel-per-pixel fitting assuming that the QSO 
emission is not spatially resolved. Left panels: residual spectra extracted 
from a region of 0.25"x0.25", where the residual map, obtained by col¬ 
lapsing the [Oiii]/15007 spectral channel, shows likely the presence of 
a spatially resolved emission. Dashed lines indicate the wavelength of 
the doublet [Oiii]/15007 and the dotted line shows the H/? position. The 
red shaded region denotes the wavelength range over which the residual 
[Oiii]/15007 emission has been integrated to produce the maps shown 
on the right. Right panels: residual maps obtained by collapsing the 
spectral channels corresponding to the residual [Om]/15007 emission, 
as shown by the red region in the left panels. In the first five maps the 
presence of a clear [Oiii]/15007 residual emission suggests that the emit¬ 
ting region is spatially resolved. The “noise” residual map of HE0251 
indicates that the sources is not resolved. 


of the [Om]/i5007 profile estimated in the preliminary spectral 
fitting. The zero-velocity wavelength was also used to refine the 
redshift of each QSO (Table[I]) providing the velocity of the host 
galaxy. The inferred redshift and [Oiii]/15007 line width are con¬ 
sistent, within the errors, with those estimated by Shemmer et al. 
(2004) and |Marziani et al.| ( |2009| ). 

The vio maps (Figure [3} third panel) show strongly blue 
shifted regions, spatially associated with high velocity disper¬ 
sion (> 400 km/s, Figure[3 forth panels). The broad [Oiii]T5007 
profile cannot be explained by a rotating gas component, which, 
in local star forming galaxies has typical FHWM values of about 
~250 km/s. Moreover, the morpholgy of the velocity maps , sug¬ 
gesting the presence of a conical blue-shifted region, is com¬ 
pletely different from the typical “spider” diagram of a disc. 
Consequently both the [Oiii]/15007 profile and the velocity maps 
suggest that in at least five out of six QSOs we detect ionised 
outflowing gas with velocities > 300 km/s. 

We can detect only the blue-side outflows because the red- 
side is likely to be obscured by dust in the host galaxy along 
the line-of-sight. For this reason, the [Oiii]T5007 line emission 
is asymmetric with a prominent blue shifted wing. While this is 
what commonly happens, in some cases the particular orienta¬ 
tion of the line-of-sight with respect to the source can result in 
redshifted outflows (e.g. Rodriguez Zaurin et al.||2013[ |Bae & 


|Woo|2014HPerna et al.|2015]T 

The next step is to understand the physical mechanisms driv¬ 
ing away the ionised gas out to a distance of a few kpc from the 
centre of the host galaxy. 


5. Outflow properties 

In the following, we estimate mass, average velocity and radius 
of the high-velocity winds. We then infer mass outflow rate, mo¬ 
mentum rate and kinetic power, and compare them with previ¬ 


ous w orks ([Greene et al.|2012[|Cicone et al.|2Q14[|Harrison et ah] 


2014, 

Sun et al. 2014, Brusa et al. 2015,|Cresci et al.P015 and 

Ferug] 

lio et al. 2015). 


5.1. Outflow and radius 


Given the uncertainties on the driving mechanism, we used a 
simple model to estimate the physical properties of the outflow. 
In this model, the outflow is represented by a shell-like cloud 
ejected from the nucleus within a cone and with a filling factor 
equal to 100%. We assume that the physical size of the outflow¬ 
ing material is smaller than the spatial resolution typical of our 
datasets. This model is defined by the mass of the cloud, M 0 , the 
distance between the cloud and the location of the AGN, R 0 , and 
the averaged cloud velocity, v Q . The outflow mass rate is given 

by 


M 0 


Mo_ 

T dyn 


M 0 v 0 


Rn 


( 1 ) 


where Tdyn is the dynamical time, i.e. the time taken by the 
ionised gas to reach a distance R 0 with an average velocity v 0 . 
We note that assuming either a shell or a uniformly filled cone 
with a filling factor equal to 1 (e.g |Maiolino et al.|2012[ ) changes 
the M 0 estimate by a factor of 3. Since this factor is constant, 
it does not affect the main conclusions of this work but we will 
take it into account when comparing our results with previous 
ones. 

The extension of an outflow (. R 0 ) is usually trivially estimated 
from the observed flux or velocity maps (e.g. |Harrison et al. 
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Fig. 3. [Om]/I5007 flux, median velocity, vi 0 , and velocity dispersion map. The maps are obtained by selecting pixels with a SNR>2. The velocity 
maps are characterised by blue shifted regions with a large velocity dispersion. Contours represent the total H/? line surface brightness at 90%, 
50%, and 30% of the peak value 
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|2014||Cresci et al.|2Q15| ). However, in our cases, the [Oiii]/15007 
emission is only marginally spatially resolved in each QSO, 
therefore the kinematical maps in Figure [3] are affected by PSF 
smearing and the sizes of the regions showing blue-shifted emis¬ 
sion (Fig. [3]) d o not directly provide the parameter R 0 to be used 
in equation [1J Furthermore, the sizes of the outflowing regions 
are affected by the choice of S/N threshold adopted in plotting 
the maps. In particular, assuming a S/N threshold larger than 2, 
the blue regions would be smaller than those shown in Figure [3] 
For these reasons, we cannot estimate R 0 from our observations 
by using the flux or velocity maps as in previous works, where 
emission lines are clearly spatially resolved. Therefore, we de¬ 
cided to estimate R 0 by using spectroastrometry, which allows 
position measurements on scales smaller than the spatial resolu¬ 
tion of the observations. Spectroastrometry consists in measur¬ 
ing the photo-centroid in each velocity channel. If outflowing 
ionised gas is moving away at a distance R from the QSO, we 
expect that the centroid of light emission extracted from blue- 
shifted velocities channels of [Oiii]T5007 will be displaced of 
the same amount R with respect to the QSO position, identified 
by the BLR and/or continuum emission. In Section [53] by means 
of a simple simulation, we show that spectroastrometry can, in¬ 
deed, provide information on the outflow position at scales that 
are significantly smaller than the limit imposed by the spatial 
resolution of the observations. 

We applied the spectroastrometry technique to the 
[Oiii]/15007 line emission after subtracting the best-fit model 
of the continuum, broad H ft and Fell emission. To maximise 
the signal-to-noise ratio for our measurements and minimize 
the uncertainty due the spectral resolution, we rebinned the 
spectra by 3 velocity channels (~ 105 km/s). The centroid of the 
[Oiii]T5007 emission in each rebinned velocity channels was, 
then, estimated by a 2-dimensional Gaussian fitting. The QSO 
position was estimated by applying the same spectroastrometry 
technique to the the continuum and the broad H ft emission. 

The results of the spectroastrometric analysis are displayed 
in Figure]?] In the left panels, we plot the distance R of the emis¬ 
sion line photocenter from the continuum one as a function of 
velocity v. The right panels show the photocenter position on 
the sky in each velocity channel. The reliable spectroastrometric 
measurements for the [Om]T5007 emission line were selected to 
satisfy the following criteria: 

o signal-to-noise larger than 1.5 for the line flux in each spec¬ 
tral channel of the rebinned spectra extracted from a 0.25" x 
0.25" 

o FWHM of the 2-dimensional Gaussian equal or larger than 
that of the PSF of the observations 

In all targets, the [Om]T5007 blue wings centroids are dis¬ 
placed at least 0.05", i.e. 0.4 kpc, from the continuum position 
and, in some cases, we observe an offset in the red wavelengths 
as well. This latter offset may be caused by three reasons: (i) 
the continuum is not well subtracted during the kinematic analy¬ 
sis; (ii) the [Oiii]/15007 emission line associated to outflows is so 
large that the displacement is slightly observable also in the red 
wings; (iii) the red wings are emitted by the receding outflows 
not completely obscured by dust in the disc. An additional expla¬ 
nation, which will be discussed in a companion paper is that in 
two out of six QSOs a fraction of the red [Ora]/I5007 line emis¬ 
sion is associated with star formation in the host galaxy causing 
a displacement from the continuum centre. 

Since the position of the photocenter at high blue velocities 
reveals the presence of the extended ionised outflow, we adopt 
R 0 equal to the largest distance measured for the approaching 
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Fig. 4. Left panels. The points show the [Om]T5007 velocity v versus 
the distance R of the [Om]45007 photocenter from the continuum (in¬ 
dicated by the dashed line). The arrows indicate the velocity v 0 = v(R 0 ) 
corresponding to the distance R 0 . Right panels. The [Ora]/I5007 pho¬ 
tocenter position in the field of view. Symbols are coloured according 
to their velocity (the velocity scale is reported in the colour bar). The 
dotted line indicates the project direction of the outflow as inferred by 
comparing the spectroastrometry results with velocity maps (Fig. [3]). 


gas, and v 0 = v(Ro) as measured from our maps (see arrows in 
Figure [4]). Since the spectra have been rebinned by 3 velocity 
channels, the error on vo due to the spectral resolution typical of 
our datasets is negligible. R 0 and v 0 for each quasar are listed in 
Table 2. 
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Table 2. Outflow properties 
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QSO 

Vo 

[km/s] 

R 0 

[kpc] 

Log(L-«fj w ) 

Log(L^ tflow ) 

M outflow 

iKZ [OIII] 

[10 7 M 0 ] 

M outflow 

m UB 

[10’ M 0 ] 

4^[oni] 
[Mg yr 1 ] 

[Mg yr 1 ] 

LBQS0109 

1850 

0.4 

43.17 

42.10 

1.2 

2.2 

60 

110 

2QZJ0028 

2300 

0.7 

43.68 

43.07 

3.8 

20 

140 

700 

HB8905 

500 

1.3 

43.95 

42.99 

7.1 

16.8 

30 

75 

HE0109 

900 

0.4 

43.75 

42.71 

4.5 

8.6 

110 

210 

HB8903 

1450 

1.9 

42.95 

41.8 

0.7 

1.2 

6 

10 

HE0251 

- 

- 

43.82 

42.55 

5.3 

6.1 

- 

- 


Notes : The outflow masses are estimated assuming a T e ~ 10 4 K and a n e ~ 500cm 3 . 


5.2. Outflow Mass 

A big challenge in estimating the mass of the wind is caused 
by its multiphase nature. Indeed, only a fraction of the mass of 
the outflows is in the warm ionised phase traced by [Oiii]/15007. 
Recent works estimate the mass of ionised outflows either us¬ 
ing Hj3 emission line (e.g |Liu et al.|2013[ [Harrison et al.|2014| ) 
or using [Om]/15007 line (e.g |Cano-Diaz et al.||2012| ). So far, it 
is not clear which is the best tracer of ionised gas powered by 
AGN feedback. To compare two different measurements carried 
out with these two different tracers, we discuss how the mass of 
ionised outflows can be constrained through the observations of 
[Om]/i5007 and H J3. 

The [Oiii]/I5007 luminosity is given by 


-[OIII] 


= f fn e n(0 2+ )j[oni](n e ,T e )dV 
Jv 


( 2 ) 


where / is the filling factor, n e the electron density, n(0 2+ ) the 
density of 0 2+ ions and j\oui\(n e , T e ) the line emissivity. n(0 2+ ) 
can be written as 


' n(0 2+ y 

\n(0) 1 

n(H) 

n(0) J 

n(H) 

n e 


n(O z+ ) = 

and, with a reasonable assumptions, n(0 2+ ) = n(0 ), 

n(0 2+ ) * (6.04 x i0- 4 10 [O/H] - [O///] °) x (1.2)" 1 x n e 

where [O/H] - [0/H] Q is the metallicity relative to solar with a 
solar oxygen abundance of [0/H] Q ~ 8.86 ( [Centeno & Socas- 


Navarro 


2008 ). The factor (1.2) 1 takes into account a 10% num¬ 


ber density of He atoms with respect to H atoms: 

n e * n(H) + 2 n(H e ) = n(H) + 2 x 0.1 x n(H) = 1.2 n(H) 

Assuming a typical temperature (T e - 10 4 K) and electron den¬ 
sity (n e - 500 cm -3 ) for the NLR the line emissivity is 

j{ oni] = 3.4 x 10“ 21 erg s _1 cnT 3 


This line emissivity was estimated making use of PyNeb ( [Lurid- 
iana et al.|2015| ). Therefore, equation ^ can be rewritten as: 


-[OIII] 


= 6.0 x io~ 4 flO [0/H] - [0/H]Q j [0m <n 2 e >V 


( 3 ) 


where < n 2 > is the volume-averaged squared density. The gas 
mass can be expressed as 


M'. 


X 


fmn(H)dV - fm p <n e > V 


( 4 ) 


where m is the average molecular weight, m p is the proton mass, 
and we have taken into count that a 10% number density of He 
atoms with respect to H atoms: 

_ m p n(H) + 4m p n(He) n{H) 

mn(H) «- n e 

n(H ) + n(He) n e 

m p (l.2)~ l n e ~ l.2m p (l.2)~ l n e - m p n e 
Finally, combining equation ^ and equation ([4]) we get: 

m p CL[ 0 iii\ 


n(H ) + 0An(H) 
n(H) + ().\n(H)' 1 '’ 


47[oiii] = 1.7 x 10 


IQ[0/H]-[0/H] q j [0///] <n e > 


( 5 ) 


where C -<n e > 2 / < n 2 >. For T e = 10 4 K and n e = 500 cm 3 
we obtain a mass of: 


M[oni] — 0.8x10 M ( 


3 ( 


C 


io [o/H]-[om 


\ I L[oui] \ I <n e > \ 1 
"/\10 44 ^rg/^/\500cm- 3 / 


Note that the gas mass is sensitive to temperature and density of 
the clouds. 

Now we similarly derive M from H (3. The H J3 luminosity can 
be expressed as 


^HJ3 


= f fn e n p jH/ 3 (n e ,T e )dV - 1.2 1 fj H p(n e ,T e ) < n 2 e > V 

Jv 


( 6 ) 


where jnp is line emissivity and n p is the proton density that can 
be written as 


n(H) 


n e = (1.2) l n e 


At the typical temperature and density of the NLR, the line emis¬ 
sivity of H J3 also estimated with PyNeb is 

j H p = 1.2 x 10 -25 erg s _1 cm“ 3 

Combining equation ([4]) with equation ([6]) of L H p we get 

m p CL H p 


Mhj3 — 0.8 

Jh (3 <n e > 

For T e ~ 10 4 K and n e ~ 500 cm -3 we obtain a mass of: 

Lh(3 \ / < n e > \~ l 
I0 44 erg/s)\500cm- 3 ) 


( 7 ) 


M H p = 1.7 x 10 9 M o C 
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Fig. 5. Cartoon showing the basic structure of our model. The outflow 
(solid blue curve) is perpendicular to the galaxy plane, v indicates the 
direction of the outflow and r is the distance from the QSO. This model 
assumes that the outflow is marginally resolved but it is not larger than 
the PSF dimension (dashed purple line). 


Comparing the two mass at the same temperature and density, 
we get 


M [oni] _ 

Mhp L H /3 


( 8 ) 


with [O/H] = [0/H] 0 , T e - 10 4 K and n e - 500 cm -3 , the same 
values which will be assumed in the following analysis. Usually 
the [Oiii]/ 15007/H/3 ratio measured in Seyfert galaxies is of the 
order of - 10 providing M\ 0 ui\/M H p - 0.5. This apparent dis¬ 
crepancy is the consequence of the different volumes from which 
[Oiii]/ 15007 and H/3 are emitted: in the above example the vol¬ 
ume of H J3 emitting gas is - 2 larger than that of [Oiii]T 5007. 
Therefore, Hj3 emission traces a larger amount of ionised gas 
mass than [Oiii]/ 15007. The mass estimated from [Oiii]/ 15007 
will then be considered as a lower limit on the ionised gas mass. 


5.3. Simulations of spectroastrometric observations of 
outflows 

In Section mo we have described our method to estimate R 0 , 
which is not affected by the spatial resolution of the data and on 
the S/N threshold used to map velocities. 

In the following we present results of simple simulations of 
outflowing gas, validating the spectroastrometric method. Note 
that the purpose of this work is not to build an outflow model 
to fit our observations, but the only purpose of this simulation 
is to show that the extension of the flux or velocity maps when 
the source is marginally resolved can provide wrong estimates of 
the outflow extension due to PSF smearing. In fact we will also 
show that an ionised clump with a given velocity and located a 
few kpc away from the AGN can appear as a blue-shifted region 
extended over -5-10 kpc in velocity maps, which is larger than 
the real outflow extension. 

We considered a simple model of an outflow (Fig. § based 
on the following assumptions. 

1. The outflow has a biconical geometry with the axis approx¬ 
imately more or less perpendicular to the host galaxy disk; 
indeed, the AGN wind, although roughly isotropic, cannot 
propagate through the galaxy disk because of its high den¬ 
sity. 

2. The receding half cone is not observed because its emission 
is absorbed by the dust in the host galaxy. This geometry best 
explains the asymmetric [Oiii]T5007 profile with a promi¬ 
nent blue wing over 1000 km/s. 

3. The surface brightness distribution of the outflow is parame- 
terised by 

„ it )model 

I(r) = 7 0 e /R ° 



Fig. 6. [OIII]/15007 emission lines extracted from a nuclear region 
of 0.25"x 0.25" . The blue line is the [Om]T5007 emission line of 
LBQS0109 after subtracting all the other best-fit components. The red 
line is the ionised emission line obtained from the simulated data. 


where the vector r is the distance from the QSO. Note that 
the surface brightness distribution and the opening angle of 
the cone are not fundamental for these simulations, since the 
spatial resolution of our data is not high enough to resolve 
the surface brightness profile. We assume an opening angle 
of 30 degrees. R™ odel identifies the photo-center of the sur¬ 
face brightness of the ionised cloud. 7 q is proportional to the 
outflow mass M™ odel defined in equation (4). 

4. The ionised gas in the outflow has an average velocity v™ odel 
and a velocity dispersion cr^ odel . 

Therefore the input parameters are the outflow mass M™ odel , the 
distance R™ odel of the ionised clouds from the centre of the QSO, 
the average velocity v™ odel and the velocity dispersion cr™ odel of 
the ionised gas. 

In order to reproduce the [Oiii]/15007 line profile of 
LBQS0109 we simulated a QSO with an outflow characterised 
by M™ odel = 1.0 x 10 7 M 0 (eq. 5), R™ odel = 0.05 arcsec (i.e ~0.4 
kpc), v™ odel = 1500 km/s and cr™ odel = 500 km/s. We also consid¬ 
ered a point source at the location of the AGN describing the gas 
at the systemic velocity with velocity dispersion - 400 km/s and 
an amplitude scaled as to reproduce the observed line profile. 
The total [Oiii]T5007 emission was then projected onto the sky 
plane, assuming a galaxy disc inclination of 10 degrees, i.e. al¬ 
most face on as typical of QSO. Finally we convolved the surface 
brightness maps in each spectral channel with the PSF and added 
Gaussian noise to match the sensitivity of our observations. 

The comparison between the [Om]/15007 line profile ex¬ 
tracted from the simulated data and that from the observa¬ 
tions is shown in Figure [6] Both spectra are extracted from a 
nuclear region of 0.25"x 0.25". The simulated emission line 
presents a prominent blue wings similarly to the real spectrum 
of LBQS0109. We then performed the kinematic analysis on the 
simulated data as described in Section pkT) The flux, velocity and 
velocity dispersion maps (Fig. [ 7 ]) obtained from the simulations 
are similar to those shown in Figure [3] extracted from observa¬ 
tions. Indeed, the velocity map shows blue shifted velocities in 
the outflow region extended over 0.5 "(i.e -4.2 kpc). So, the sim¬ 
ulation confirms that the blue-shifted region suffers from beam 
smearing and its size does not match the real radius of the out¬ 
flow. 

Since the simulated outflow well describes the asymmetric 
[Oiii]/15007 profile and the velocity gradient observed in our 
datasets, we carried out the spectroastrometry analysis on the 
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Fig. 7. [OIII] flux, velocity and velocity dispersion map obtained from 
the simulated [Om] 45007 data. The regions with S/N < 2 are masked 
out. 


simulated data to estimate the average velocity and radius of 
the ionised gas in order to compare these derived values with 
the input parameters. The offset of the photo-centroids at differ¬ 
ent velocities are shown in Figure [8] The radius R a (~ 0.05 ") 
and the average velocity v(R 0 ) (~ 1450 km/s) estimated from 
spectroastrometry are consistent to the initial input parameters. 
So, provided that the S/N is high enough, the spectroastrom¬ 
etry method can provide the position and the velocity of out¬ 
flows with an accuracy well below the seeing limit. Moreover, 
since the integral flux of [Oiii]/15007 at high blue-shifted veloc¬ 
ity (-2500km/s < v < v(R 0 )) is consistent with about half of the 
total flux emitted from the ionised outflowing gas, the outflow 
luminosity can be inferred from: 

L 0 = 2 [°L A (A)dA (9) 

where To is the wavelength corresponding to vo. In the cases 
where the [Oiii]/15007 emission line was well described by a 
multi-Gaussian fit, the luminosity value calculated with equation 
0 and the one estimated from the broad Gaussian component, 
as typically done in the literature, are consistent within the er¬ 
rors. So, in the multi-Gaussian fit we estimated the L 0 from the 
broad component and in the other case we measured the lumi¬ 
nosity of the gas using the equation 0- 

Since the narrow H (3 component has the same profile of 
[OIII], we can measure the outflow luminosity for both the emis¬ 
sion lines using the same equation. From the luminosity we can 
get the mass of the ionised outflows by using equation and 
([7]). As the [Oiii]T5007/H/3 ratio measured in each QSO is ~ 10, 
the mass inferred by Hj3 is ~ 2 times larger than that estimated 
by [Oiii]T5007 (see eq. [8b. The luminosity and mass of the out¬ 
flow are listed in Table jl] and [2| respectively. Since HE0251 is 
not spatially resolved, we cannot estimate the outflow luminos¬ 
ity using the method described above. In this case we fitted the 
[Oiii]/15007 with two Gaussian components and we inferred the 
outflow luminosity from the broad Gaussian component. 


6. Results 

In order to investigate the nature of the ionised wind, we need 
to estimate the main quantities of the outflows: mass outflow 
rate, momentum rate and kinetic power. Theoretical models (e.g 


Zubovas & King|20T2] |Faucher-Giguere & Quataert|2012| ) pre¬ 

dict tight relations between these quantities and the AGN bolo- 
metric luminosity, Lagn- In particular the correlation between 
the momentum rate and L A gn provides an indicator of the na¬ 
ture of the feedback mechanism. In this section we derive the 
main properties of the AGNs in our sample and compare them 


with the prediction of the models. As explained in Section 5.3 


the low spatial resolution of our observations does not allow us 
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Fig. 8. Spectroastrometry results obtained from the simulated data. On 
the left, centroid offset of [Oiii]T5007 emission respect to continuum 
emission (dashed line) at different wavelength. On the right, the photo¬ 
centroid position in each velocity. Colour coding corresponds to veloc¬ 
ity offset. 


to discern whether our outflows are single explosive events or 
refilled with clouds ejected from the galactic disk. Hence, we 
assume a simple model where a single ionised cloud is ejected 
outward of the nuclear region and the mass outflow rate is equal 
to the mass of the outflow divided by the dynamical timescale 
(see eq. 1). The dynamic time is the time that a clump of ionised 
gas in outflow takes to reach a fixed distance from the QSO. We 
derive M 0 values in the range 6- 700 M© yr _1 . We assume an 
average outflow mass rate error of ±50%, which takes in account 
the uncertainties associated with the gas physical properties (i.e. 
density and temperature), flux calibrations and Ro and Vo esti¬ 
mates due to projection effects. For each source we estimate two 
value of M 0 that are calculated by using Hj3 and [Oiii]/15007 re¬ 
spectively. Since the [Oiii]/15007 outflow luminosities are ~ 10 
times larger than the ones from Hj3 , the outflow masses inferred 
from the hydrogen emission line are ~2 times larger than those 
estimated by [Oiii]/15007. 

In Figure [9] we plot the outflow velocity and mass rate as a 
function of the AGN luminosity. Lagn is derived by using the 
relation L A gn ~ 6TL(T5100A) from Marconi et al. (2004). The 


solid and hollow blue points are the M 0 estimated for five out of 
six QSOs; we could not infer the value of HE0251 since we were 
not able to measure the size of the outflow. The red points rep¬ 
resent the ionised outflows observed in type 2 AGN at redshift 


0.08 < z < 0.2 (Harrison et al. 2014 

); the orange and purple 

points are from 
mapped ionisec 

Brusa et al. ( 2015) and 

Cresci et al. (2015 ) who 

outflows in 6 X-ray selected, obscured QSOs 


at z 


1.5; the black circle corresponds to the the ionised out¬ 
flow in a obscured radio-quiet QSO at z -0.123 ( Gree ne et al. 


2012). For consistency with our work we re-calculated the out¬ 
flows properties using an electronic density of n e =500 cm -3 and 
a temperature of T e = 10 4 K. Note that the outflow masses from 
literature may be overestimated since we used the luminosities 
from the total [Oiii]T5007 and H (3 profiles. 

In addition to ionised outflows, we compare our results with 
those obtained from molecular outflows (green, black and brown 
solid squar e; |Cicone et al.||2014| |Sun et al.' |2014| and |Feruglio| 
et al.|2015| respectively). The molecular outflow properties were 


re-estimated assuming a shell-like cloud model (eq.[T]). 

In Figure [9] we show that the outflow velocity vo and out¬ 
flow rate are correlated with the AGN luminosity, although with 
a large scatter. The increase in velocity and outflow rate with 
increasing AGN luminosity is consistent with the idea that a lu¬ 
minous AGN pushes away the surrounding gas through a radia- 
tively driven fast wind whose kinetic power is a fraction of the 
AGN luminosity. However, it is not easy to establish a direct re- 
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Fig. 9. Left: Outflow velocity as a function of the AGN bolometric luminosity. The blue circles denote the results from this work; the open circles 
mark the velocities of ionised gas (mainly [Oiii]/15007): the red, orange, purple and black circles are the estimates obtained from Harriso n~et al.| 
p014| ), [Brusa et a l. ('2015), Cre sci et al.|(|2015|) and|Greene et al. ( 2012), res pectively. The green , black and brown squares denote the velocities 
of the molecular outflows from |Cicone et al.| ( |2014| ), |Sun et al. 1(2014) and fFeruglio et al.| < [2015]) . Right: Outflow rates as a function of the AGN 
bolometric luminosity. Notation is the same of the left panel expect that open circles mark the estimates obtained with [OIII]-inferred masses, 
while the filled ones denote the estimates based on H fi. We recalculated the outflow properties inferred by previous works to make them consistent 
with our estimates (see text). The solid line is the best fit relation to the averages of filled and empty circles. 


lation between AGN luminosity, velocity and outflow rate. The 
acceleration process of the outflow by the fast wind and the frac¬ 
tion of the kinetic power injected, might vary from object to ob¬ 
ject; moreover, the observed AGN luminosity may not represent 
the long-term average luminosity which is at the end respon¬ 
sible for driving the outflow. For these reasons we do not ex¬ 
pect tight correlations between velocity, outflow rate and AGN 
luminosity, as observed. Indeed, previous studies on more het¬ 
erogenous samples, but smaller luminosity ranges, did not find 
any significant correlation between the outflow velocities and the 
AGN luminosities (e.g. |Veilleux et al.|2013||Brusa et al.|2015| ). 
In the luminosity range where we have data from both CO and 
[OIII] (10 44 - 5 erg/s < L A gn ^ 10 46,5 erg/s), the velocities of the 
ionised outflows are a factor ~2 larger than those of the molec¬ 
ular gas, consistently with larger masses of molecular gas in 
both momentum- and energy-driven scenarios, but also with dif¬ 
ferent acceleration/deceleration processes. However, given the 
heterogeneous nature of the sample and the non-uniform mea¬ 
surements of outflow velocities we cannot draw any firm con¬ 
clusions. Moreover, there is only one galaxy (black circle and 
square in Fig. [9]) where both molecular and ionised outflows are 
detected. This discrepancy is even more evident when comparing 
molecular with ionised outflow rates. |Cicone et al.| ( |2014| fit a a 
log-linear relation between M and L A gn for their molecular out¬ 
flows finding Logio(M) = 2.84 + 0.720 xL<9£io(Lagv/10 46 'erg/s) 
where we have corrected for the factor 3 discrepancy described 
in Section 5.1 Here, we fit the same relation between M and 


Lagn fixing the slope to that of |Cicone et al.| ( |2QT4] ): 


Logio(M) = (0.6 ± 0.2) + 0.720 x Lo^i 0 (L AGiV /10 46 erg/s) (10) 


Since the velocities and the radii of ionised outflows are simi¬ 
lar to those of molecular ones, we interpret the offset between 
the two Lagn- M correlations as an indication that the ionised 
gas only traces a small fraction of the total gas mass. Comparing 
the normalization of the Lagn- M relation inferred from our data 
with that obtained from molecular outflows, we can infer that 
the outflow rate of the ionised gas is a small fraction (< 10%) 
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Fig. 10. Kinetic power as a function of the AGN bolometric luminos¬ 
ity. Symbols and colours as in fig. [9] The solid, dashed and dotted line 
correspond to P k = 100%, 5%, 0.1 %L AGN respectively. 


of the molecular one, suggesting that outflow masses are dom¬ 
inated by molecular gas. However, since molecular and ionised 
outflow rates are not estimated in the same objects and the rela¬ 
tions have quite a large scatter, it is not possible to estimate the 
exact fraction of ionised gas in the outflows. On the other hand, 
the different normalisation could indicate that two different ac¬ 
celeration mechanisms are at work. Future ALMA observations 
of the molecular gas in these objects are essential to distinguish 
between these two different scenarios. 

Figure [10] shows the kinetic power of outflowing gas as a 
function of AGN luminosity where the outflow kinetic energy 
rate is given by 
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Fig. 11. Outflow momentum rate as a function of photon momentum of 
the AGN . Symbols and colours as in Fig. [9] The dotted and dashed line 
correspond to 20L A gn/ c and P~L AG n/ c respectively. 
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We indicate using solid, dashed and dotted lines the locus of 
points having an outflow kinetic power that is 100%, 5% and 
0.1% of the AGN luminosity, respectively. The recent AGN feed¬ 
back models (e.g. |King|2010[ |Zubovas & King 2012; Lapi et al.] 
2014 ) predict a coupling efficiency between AGN-driven out¬ 
flows and AGN power of about ~5% which is needed to explain 
the M bh - cr relation observed in local galaxies. Although the 
molecular outflow observations are consistent with the models 
within the error, the kinematic power estimated from ionised out¬ 
flows is only < 0.1% of the L A gn- This percentage is too low to 
explain the M B h - cr relation. Since the kinetic energy is propor¬ 
tional to the outflow mass, we cannot discern whether the nature 
of ionised outflows is different with respect to molecular ones 
or, as discussed above, the [Om]/15007 (and H J3) line emission 
traces only a small fraction (< 10%) of the total outflowing gas. 

The last fundamental parameter of the outflows is the outflow 
momentum rate defined as v 0 M 0 . Energy-driven outflow models 
predict that the momentum of the large-scale outflow is boosted 
compared to the nuclear wind (or AGN radiation pressure mo- 
mentum ), i.e. v 0 M 0 ~ 20L A gn/ c (Zub ovas & K ing 2012). In 


Figure[TT]we observe that molecular outflows 
predicted by models (see Cicone et al. 2014 


ollow the relation 
for more details) 
while ionised outflows are below the 1:1 relation. This could 
simply be the consequence of the discrepancy between the out¬ 
flow masses of ionised and molecular gas. 

Finally, we may conclude that the main difference between 
molecular and ionised outflows is the fraction of total gas mass 
pushed away from the AGN-driven wind. This is confirmed by 
the larger outflow rates (right panel of Fig. [9]) and kinetic power 
(Fig. [TO] ) of molecular outflows, and by their smaller outflow 
velocities (left panel of Fig. [9]). 


7. Conclusions 

We performed seeing-limited, near-IR integral-field spectro¬ 
scopic observations with SINFONI of a sample of six, high lu¬ 
minosity (Lboi = 10 13 - 10 14 L 0 ) QSOs at redshift z ~ 2.4. 
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[Oiii]/15007 emission lines, redshifted into the H-band, are char¬ 
acterised by large FWHMs (> 1000 km/s) and prominent blue 
wings, indicative of fast outflows accelerated by the powerful 
AGN. We find that: 

o The [Oiii]T5007 emission line is spatially resolved in five out 
of six sources and extended over several kiloparsecs. 
o The analysis of the [Oiii]/15007 kinematical maps suggests 
the presence of conical outflows associated with regions of 
high velocity dispersions ( > 500 km/s) The physical prop¬ 
erties of the outflows, i.e. mass, outflow rate, kinetic energy 
and momentum rate, have been estimated with a new method 
based on spectroastrometry which is not affected by PSF 
smearing, at variance with results from previous works. The 
reliability of this method has been confirmed with a simple 
simulation but more accurate modelling will be developed in 
future work. 

o Both ionised and molecular outflow velocities are weakly 
correlated with the observed AGN luminosity, despite the 
large scatter, and at a given AGN luminosity the velocity of 
the ionised gas is roughly a factor 2 larger than that of the 
molecular one. However, given the heterogenous nature of 
the combined sample we cannot draw any firm conclusions, 
o Mass and momentum rates, as well as kinetic powers, in¬ 
crease with AGN bolometric luminosity in a similar way 
to what is observed in molecular outflows in the local uni¬ 
verse. The ionised gas properties define relations with AGN 
luminosity which are parallel with those of molecular gas. In 
particular, ionised outflow rates are ~ 50 times lower than 
molecular ones. The kinetic power carried by ionised out¬ 
flows is of the order of ~ 0.1 - 0.05% of the AGN luminos¬ 
ity compared to ~ 5% for molecular outflows. Finally, mo¬ 
mentum rates are of the order of L AG n/c, a factor ~ 20-50 
smaller than for molecular outflows, 
o These discrepancies between ionised and molecular outflows 
can be explained with the fact that ionised gas traces a 
smaller fraction of the total gas mass. Alternatively, they 
are the indication of different acceleration mechanisms for 
the molecular and the ionised gas. Observations with AFMA 
would allow us to measure the molecular gas mass in these 
objects therefore distinguishing between these two scenarios. 
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